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Summary 

A large literature exists concerning the equilib- 
rium configurations of free liquid/gas surfaces 
in reduced gravity environments. Such condi- 
tions generally yield surfaces of constant curva- 
ture meeting the container wall at a particular 
(contact) angle. The time required to reach and 
stabilize about this configuration is less studied 
for the case of sudden changes in gravity level, 
e.g. from normal- to low-gravity, as can occur 
in many drop tower experiments. The particular 
interest here was to determine the total reorienta- 
tion time for such surfaces in cylinders (mainly), 
as a function primarily of contact angle and kine- 
matic viscosity, in order to aid in the development 
of drop tower experiment design. A large para- 
metric range of tests were performed and, based 
on an accompanying scale analysis, the complete 
data set was correlated. The results of other in- 
vestigations are included for comparison. 

X. Motivation and Background 

The use of drop towers to study liquid behavior 
in low gravity environments is common. Due to 
the step reduction in gravity level, a character- 
istic of most towers, reorientation of free liquid 
surfaces often occurs upon entry into low grav- 
ity. For cylinders, the degree of reorientation is 
dependent on the contact angle (0) of the partic- 
ular liquid/solid pair and the change in system 


Bond number (J9o), where 


with p , a, R y and a denoting liquid density, ac- 
celeration due to gravity or other force, cylinder 
radius, and fluid interfacial tension, respectively. 

In general, the Bond number of equation (1) is 
a ratio of body forces to capillary forces and can 
change several orders of magnitude during the 
drop release process. In cases where the forces 
influencing free liquid interface shapes shift from 
the gravity dominated case (Bo > 1) to the situ- 
ation where capillary forces dominate (Bo < 1), 
surface reorientation is most pronounced upon 
entry into low gravity; the fluid will flow from a 
“flat” configuration towards its low gravity equi- 
librium shape, one generally of constant curva- 
ture meeting the container wall at the contact 
angle of the system as pictured in Figure 1. De- 
pendent on the damping of the system, an over- 
shoot can result and oscillation about the equi- 
librium surface shape occurs until damped. In 
many cases the normal and low gravity equilib- 
rium surface configurations are determinable [1]. 
However, the transients associated with the re- 
orientation have not as of yet been quantified. 

A number of related analyses have been per- 
formed concerning surface reorientation and/or 
oscillation in low gravity environments [2-13]. A 
main thrust of these studies has focused on exter- 
nally forced systems and the problem of lateral 
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sloshing. In addition, fluids and tank geometries 
were selected which allowed legitimate inviscid 
flow approximations with a nearly 0° contact an- 
gle condition at the wall. Most specifically, refer- 
ences [9] and [12] approach the objectives of this 
study. Unfortunately, the former used low vis- 
cosity fluids and did not have ample low gravity 
time to determine a total settling time. The lat- 
ter used a similar experimental approach as that 
of this study, but to observe a distinctly unre- 
lated phenomenon. In a Bmall set of preliminary 
experiments Ross and Pline [13] studied axial os- 
cillations and determined qualitatively that sur- 
face oscillations and damping were strong func- 
tions of viscosity and contact angle. This effort 
expands the data base of these earlier studies, 
quantifies the settling times for a number of liq- 
uid/container (with air) systems, and presents 
functional relationships based on a scale analysis 
to determine the general settling times for other 
cylindrical systems. 

2. Scale Analysis 

A scale analysis was performed to identify the 
important parameters for which the surface set- 
tling time data of this study could be presented. 
The general form of the scaled parametric equa- 
tions was then used for correlation purposes. 

The notation and language of this section fol- 
low those of Bejan [14], where commas between 
parameters denote some functional relationship, 
normally addition, while the similar sign (~) 
denotes “of general order”. Actual operators 
(x ,+...) are retained only when known. Divid- 
ing by scaled quantities often demands physical 
insight into the problem since, for example, the 
scale function 1 /(a 2 , 1) could also be interpreted 
as (a _2 ,l). Scale constants appearing in equa- 
tions are considered 0(1) unless they arise from 
some known relationship. 

2.1 Scaled Momentum Equation. 

The analysis begins with the momentum equa- 
tion written in cylindrical coordinates for the ax- 
ial component, z, 


d lL + V d lL + V d A i = + 

dt + r dr + ‘ dz pdz 


(d 2 V, IdV. PVA 
\ dr 2 r dr dz 2 ) 


+ ng 


( 2 ) 


where ng represents the body force term, n be- 
ing the ratio of the low gravity level to normal 
earth gravity 1 . Direct scaling of each term in 
the above equation follows when an axial length 
scale is determined. This scale, Z } is selected as 
the axial displacement of the meniscus between 
an assumed initially flat interface and a steady 
end state surface condition of constant curvature 
(n = 0). Figure 1 depicts Z geometrically as a 
function of the contact angle (0) and the radius 
( R ) of the cylinder. 

A time Beale, tf f for the momentum equa- 
tion is defined as a “flow time” or characteris- 
tic response time for the system. By introduc- 
ing the scale relations V x ~ and —dP/dz ~ 
a cos 9/RZ y and by eliminating Z for the geomet- 
rically equivalent aiZ, the momentum equation 
may be rewritten as the scale equation 


aR c cos 0 

Tf* ~ paR 2 ’ 


v 

atfR 


(a 2 , 1) , ng 


where 


1 — sin 0 
cos 0 


( 3 ) 

( 4 ) 


and ranges 0 < |a| < 1. Grouping the time- 
dependent terms of equation (3), the requirement 
for which ng may be ignored is the case of Bo n <C 
1, or 


Bo n 


pngR 2 
a (1 + sin0) 


( 5 ) 


where Bo n is a contact angle-dependent, low 
gravity Bond number. This condition is satis- 
fied for all of the tests performed in this study 
and for a great many other drop tower experi- 
ments due to a general size limitation. As a re- 
sult, the body force term need not be considered 

1 Neglecting the body force term, os will be accom- 
plished shortly, produces similar scale results os would 
the anolysis if it had begun with the momentum equation 
for the rodial component. Axisymmetric reorientation of 
the fluid is ossumed, or V = T(r, z t t). Also, when lim- 
ited to the liquid/gas interface scenorio, the momentum 
balance for the gaseous phose moy be ignored. 
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in the scale analysis except to determine the ini- 
tial condition of the surface. Thus, the scaled 
terms of importance are the unsteady, capillary, 
and viscous terms, or 


aR a cos 6 v 


tj 1 pR 2 a 9 if Rot 


(i.« 2 ) 


( 6 ) 


respectively. The time scale, £/, may then be de- 
termined for the particular cases as listed below: 

unsteady ~ capillary ; when ( (l, a 2 ) < 1 


t tl ~ — Co 2 (?) 

v 

capillary ~ viscous ; when ( 2 (l,a 2 ) 1 

o') (8) 

unsteady ~ viscous ; itfhen ( (l,a 2 ) 1 


associated with “nearby” walls, which, by sig- 
nificantly altering the shape of the surface, can 
also greatly influence the dynamic behavior of the 
fluid. 

2.2 Scaled Energy Equation. 

The desired settling time scale, t B} will be some 
function oft/. To determine this relationship, a 
general end state energy balance may be written 

AE g + AE* = (13) 

where subscripts g and <j denote the change in 
hydrostatic and interfacial potential energies re- 
spectively. Efi represents the viscous dissipation 
by which the potential energy is converted to 
heat. Accurate theoretical values for the two left 
hand quantities may be determined. However, 
E M is inherently a scaled term in this analysis, 
therefore, each term is presented in a generalized 
scale form: 

AE g ~ (, 0pgR ) 2 H 3 fp, \pgR 2 H 2 (14) 


if. 


R 2 a 2 

V (1,Q 2 ) 


where 



t ) 

Rcra 2 cos 0/ 


1/2 


( 9 ) 


AE„ ~ 2nR 2 ficr cos 0 (15) 


( 10 ) 


l*dVdi~^f 2 f a j‘'-± jdt (16) 

where 


When the condition of £ 1 is met, the expres- 

sions of equations (8) and (9) are valid for tf. 
Therefore, 




V 


R 2 a 2 

v (l,a 2 ) 


( 11 ) 


fp=P 2 +P+\ (17) 



Ai 

27t R 2 cos 6 


(18) 


or 

~ ^ t<’ (*. <*=) . 11 <«) 

Since ( > 1, ( 2 » 1 and the contribution of t/ # 
in equation (12) is negligible, and subsequently 
ignored. This leaves only tf x and for consid- 
eration. 

The general forms of equations (7) and (8) are 
well known [15]. The uniqueness of these results, 
however, is the inclusion of contact angle effects 



jfe, /?, and p represent the fluid elasticity (units 
of pressure), compressibility coefficient (V _1 5 V 
/dP) t and dimensionless reference pressure, 
P 0 /pgH , respectively [16], with P a being the am- 
bient pressure of the gas phase. H is the height 
of the liquid column and the volume of the fluid 
is denoted by V. The functions f\ and /2 result 
from geometric considerations. Specifically, f\ 
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arises from the known change in interfacial area 2 
and /2 from the volume of fluid involved in the re- 
orientation. Ai is the initial liquid/gas interfacial 
area and is a function of 9 and the initial, nor- 
mal gravity Bond number (Bo* = pgR 2 / a). f a is 
a collection of the scaled nonlinear terms of the 
dissipation function ($) and may be expressed as 

fa ~ (1, a 2 > a 4 ) (20) 


A Eg contains energies associated with fluid com- 
pressibility and the initial normal gravity pres- 
sure gradient of the standing liquid column. The 
former may be neglected since the fluids of this 
experiment are virtually incompressible (/? — ♦ 0). 
Comparing the relative importance of A E g (with 
(3 = 0) to A the Bond number, Bo#, based 
on liquid column height is defined as 


Boh = 


A 

A B 9 


pgH 2 

ficr cos 9 


( 21 ) 


For all cylinders we studied, Boh 1, imply- 
ing AE g A However, the experimentally 
observed settling time was independent of H . 
Subsequent flow visualization via particle seeding 
and a laser light sheet revealed that the fluid mo- 
tion was limited to the near-surface region con- 
tained within the height Z. These two observa- 
tions suggest that the potential energy stored as 
gravitational potential (AE g ) is not dissipated by 
viscosity, but by some other mechanism 3 . Subse- 
quently we remove A E g from the scale equation 
(13) which yields 

A B. ~ (22) 


2tT j R 2 /i<7COS0 ^ 



(23) 


2,3 Settling Time Scale, t,. 

Substituting the flow time scales of equations 
(7) and (8) into (23) separately, and noting 

2 This area Include* all phases, both wetted and non- 
wetted surfaces. 

3 This explanation is in direct contradiction to that of 

ref. [12]. 


f*'dt = t t) two Bcaled settling time relations may 
be expressed with the provided constraints 


t 




R 2 ficc 4 

*hfa ’ 


<(l,a 2 ) < 1 


(24) 


(l.o, 3 )’ > 1 (2B) 

V J2J a 

The constraints on equation (24) and (25) allow 
for the superposition of time scales t it and 
leaving 

~ l » 2 


= — 77- (1, «’)■’] < 2 «) 

^ J 2 Ja L J 

When Boi > 1, implying an initially flat liq- 
uid/gas interface (j4x = *R 2 ), the quantity f\j f 2 
reduces to a constant 0(1). This final assump- 
tion noted, equation (26) reduces to 

(l.« 2 )' 2 ] (27) 

Vja. 1 J 

which may be more clearly presented as 

t. ~ ^ f(0 ) _1 [C 2 . 9(9)] (28) 

or nondimensionally as 

r= t gm~m, gw ( 29 ) 

This suggests that of the five variables of the 

problem established (R t v, p, < 7 , 0), only two pa- 
rameters, ( and 0, need be varied experimentally 
to provide ample data to correlate r. This is 
accomplished primarily by varying the fluid vis- 
cosity and contact angle in the experiments to be 
discussed shortly 4 . 


4 The important dimensionless groups may be identi- 
fied in their more common form by dividing eq. (27) by 
multiplying through by and disregarding con- 

tact angle dependence, yielding 

Rvp R 2 

1 'V — , 

<rt, ut t 


Introducing the scale relation V ~ R/t» and substituting 
into the above yields 


fjV VR 

<T V 


~ Co, fie 


where Co and He are Capillary and Reynolds number 
equivalents. This dimensionless number dependence is 
anticipated [17]. 
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The complexities of f(6) and j(0) in equation 
(29) need to be addressed with insight gathered 
from experiment, since, for example, the scale 
form of /(<?) 


m~ 



(l,° 2 ,a 4 ) 

q 4 ( 1 , a 2 ) 2 


(30) 


may have numerous interpretations. 


tube/syringe technique was adopted to reduce 
this effect as is shown in Figure 2. Just prior 
to the drop a prescribed amount of fluid was in- 
jected and the syringe was removed. The lights 
and camera were initiated approximately 1.5 sec- 
onds before the drop. The release of the package 
began a delay /relay operation which maintained 
power to the lights and camera during the drop. 


3. Experiment 

3.1 Apparatus. 

A schematic of the experimental apparatus is 
shown in Figure 2. The main components are 
a circular cylinder whose axis is oriented with 
the gravity vector, a leveling pad, top and back 
lighting, and amotion picture camera. The cylin- 
der diameter used for the bulk of the testing was 
19.05mm to ensure complete settling within the 
drop time (less than 2.2 seconds). 

The cylinders were drilled out of a rectangu- 
lar bar of acrylic plastic to minimize refraction 
associated with the vessel curvature and the mis- 
matched fluid/acrylic refractive indices. All plas- 
tic surfaces were polished to optical clarity with 
progressively finer grit SiC paper and a buffing 
compound. 

For most tests the internal surfaces of the cylin- 
ders were prepared by washing them twice in 
refrigerant-12. They were then rinsed twice with 
ethanol and cleaned ultrasonically in a warm mild 
detergent for 30 minutes. After the cleaning, the 
cylinders were again rinsed with ethanol and al- 
lowed to air dry. In this way the surface condi- 
tions of the vessels and resultant data were made 
repeatable. 

Poorly leveled cylinders produced lateral slosh- 
ing modes easily observed on the films. To reduce 
this effect a leveling pad was employed. The pad, 
made of 1cm foam sheet, was squeezed between 
the cylinder base plate and package floor. Sus- 
pended from a wire, the entire experimental ap- 
paratus was balanced and the cylinder righted via 
iterative adjustment of the wingnut fasteners. 

Another cause of slosh/swirl during the ex- 
periments was attributed to an asymmetric pre- 
wetting of the cylinder walls which occurred if 
the fluid was disturbed (sloshed) before the re- 
lease of the package in the tower. A simple fill 


3.2 Fluids Selection and Experiments. 

The 200 Fluid series (varying viscosity silicone 
oils) produced by the Dow Corning Corporation 
and aqueous glycerine solutions allowed variation 
of the dynamic viscosity without significant vari- 
ation in other relevant fluid properties. In order 
to vary only the contact angle with the silicone 
oils retained as test fluids, a thin film of surface 
modifier (F-723 supplied by the 3M Co.) was 
applied to the cylinder walls. The film was ap- 
plied just prior to each drop and just after each 
cleaning procedure by filling the cylinder with the 
surface modifier as a liquid and slowly draining 
it from the bottom. Two coats of F-723 were 
applied to each cylinder for these tests and the 
“coated” cylinder was allowed to air dry. The 
relevant properties for the silicone oils and other 
fluids are listed in Table 1. As seen from the table 
three decades of £ (with 0 = 0°) were obtained 
by use of 0.65 through 200.0cSt oil. 

A series of tests was also performed using 
aqueous ethanol mixtures in un-coated cylinders. 
Since distilled water and ethanol produce equi- 
librium contact angles of 80° and 0° on acrylic, 
respectively, combinations of the two fluids were 
selected to produce a range of angles. The sur- 
face tensions, viscosities, and densities of these 
solutions and those of the aqueous glycerine solu- 
tions were measured specifically for this research 
work using standard laboratory procedures and 
are shown in Table 1 with their respective uncer- 
tainties. The equilibrium contact angles for each 
fluid/solid pair were measured using a number of 
techniques: the capillary tube, a modified Wil- 
helmy technique, and a sessile drop (tilt-slide) 
method. Because contact angle hysteresis is a 
problem for each of these methods, an equilib- 
rium contact angle was determined by a fourth 
technique exploiting the steady, end state config- 


5 



uration of the low-g liquid surface as suggested in 
[18]. Since the surface forms a shape of constant 
curvature, a geometric relationship was derived 
relating the axial length scale, Z y to the contact 
angle (ref. Fig. 1), 


a = — 

R 


1 — sin 9 
cos 9 


(31) 


The contact angle may then be expressed as 


9 = ,in '' (irfs) 


(32) 


where a is determined experimentally. Equilib- 
rium contact angles determined in this manner 
are listed in Table 2. In each case the equilib- 
rium angles were found to be within the range 
of hysteresis determined by the above mentioned 
techniques. 

Supplementary tests were performed to note 
geometric effects and the dynamic behavior of the 
contact line during the reorientation and settling 
of the surface. Concerning the former, square 
vessels were employed for these tests, and, due 
to their larger characteristic length, the 5 sec- 
ond drop facility at NASA Lewis Research Cen- 
ter was used. Fluid selection was made similarly 
as with the cylinders in order to vary primarily 
the fluid viscosity and contact angle. As is the 
case for capillary flow in corners, stable surfaces 
result when the condition 9 > 90° — ip i s satis- 
fied [19], where ip is the half-angle of the corner in 
question. For the Square vessel ip = 45° and con- 
tact angles in this range were tested. As a point 
in verification of this criteria, contact angle val- 
ues which did not satisfy this criteria produced 
surfaces which, after having “climbed” or wicked 
up the corners, formed relatively complex shapes 
or wrapped entirely around the ullage. Unfor- 
tunately, the results of these tests were predomi- 
nantly qualitative as positive identification of the 
interface was always in question due to a light- 
ing and photographic limitation. It can be said, 
however, that the general settling behavior for the 
vessel with square cross section followed similarly 
that of the right circular cylinder. 

The transients of the settling process could 
not be observed at the contact line for the ma- 
jority of the tests. This was attributable to 


the significantly mismatched refractive indices of 
the test fluids (n = 1.35) and acrylic contain- 
ers (n = 1.49) employed. In order to view the 
contact line region, a special fluid solution was 
prepared possessing an identical index of refrac- 
tion as that of the acrylic. The solution consisted 
of Cargille Laser Fluid (n = 1.5) and 2c5t sili- 
cone oil (n = 1.39) and was iteratively mixed and 
measured using a standard refractometer until it 
matched the value of the acrylic. The precision of 
the measurement under controlled temperature 
conditions was better than 0.5%. The viscos- 
ity, density, surface tension, and contact angle on 
an FC-723 coated surface were also measured for 
the mixture and found to be 14.9c5t, 880Ay/m 3 , 
26.3AT/m, an <3 70°, respectively. The dependence 
of the solution refractive index on temperature 
was also determined and, for the range of tem- 
peratures at which the drop tests were performed, 
was shown to not produce significant changes in 
the value. 

4. Results and Discussion 

4.1 Primitive Variable Dependencies. 

Selected surface settling time histories of the 
meniscus apex for a given cylinder diameter are 
plotted in Figures 3 and 4 for the silicone oil tests. 
Figure 3 is presented to display viscous effects 
( 9 — constant) while Figure 4 is to show con- 
tact angle effects ( u = constant). The data for 
these figures were obtained from the film records 
by taking measurements, for the most part, ev- 
ery 1/64 sec to ±0.05mm. Data read at 1/128 
sec produced identical results. Settling time and 
frequency characteristics were directly obtained 
from such data and are compiled in Table 2 for 
the various tests performed. 

The frequencies of the surface oscillations were 
found to be regularly periodic and were de- 
termined by averaging peak-to-peak values-the 
number of peaks (nodes) used in the averages are 
listed in Table 2. The settling time was defined 
as the point at which regular oscillations were 
no longer detectable or, in the case of the high 
viscosity fluids experiencing creep, the point at 
which the axial meniscus location became steady 
in time. In some cases these settling characteris- 
tics were additive. This is exemplified in Figure 
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3 for the lO.OcSt silicone oil test. The settling 
time for such a case was defined as 1.75sec and 
not 0.75sec. 

Due to the high magnification of the camera on 
the cylinder, the fore-field of view (leading edge 
of the meniscus) was blurred yielding a “thick” 
initial interface 0.3mm) which was difficult to 
identify consistently. As a result, the initial sur- 
face position varied to the extent that the equilib- 
rium surface location did not always coincide for 
identical contact angles (refer Figure 3). How- 
ever, frequency and settling time measurements 
were not affected by this condition. 

The dynamics at the contact line during the 
process are illustrated in Figure 5 for the index- 
matched solution. It is evident from the figure 
that oscillation at the contact line is not present. 
It should be stated, however, that due to the re- 
fractive index limitation, fluid viscosities below 
xcSt could not be tested and non-oscillatory be- 
havior of the contact line, or even contact angle, 
should not be assumed for all values of (. For 
the case of Figure 5, variance in the contact an- 
gle could not be determined. 

The raw frequency and settling time data for 
the silicone oil tests is presented in Figures 6 and 
7 for a cylinder of radius 9.52mm. The damping 
ratio with 8 = 0°, or £ 0 = 0 ° = v(p/Rc r) 1 ^ 2 , is 
used on the abscissa of both figures. This is done 
to display the effects of varying contact angle on 
the resulting settling characteristics. It should 
be noted that the contact angles for the silicone 
oil/F-723 tests ranged from 32° to 67° (points 
left to right on plots in order as listed in Table 
2). Regardless of this fact, a strong contact angle 
effect is noticeable. The settling time data for the 
silicone oil/F-723 tests suggests an upward trend 
for large values of £ as does that of the silicone 
oil/acrylic system (0 = 0°) at ( = 0.003. 

4.2 Nondimensionalized Experimental Results. 

The data of Figure 7 is again displayed on 
Figure 8 with t, nondimensionalized and plot- 
ted against ( following the form of equation (29). 
The function f(8) = 1 + a 2 is used so that 

r = % (> + « ! ) (33) 

Also shown on Figure 8 is the aqueous ethanol 


data set with R = 9.52mm. The primary vari- 
able for the aqueous ethanol mixtures was the 
contact angle which varied from 5° to 67° , the 
lower values collapsing on the 8 = 0° silicone 
oil/acrylic curve. From this figure three char- 
acteristics of the r-( relationship were noted for 
the range of ( tested: 1) for 8 held constant, r 

obeys a power law dependence on 2) for con- 
stant £, r decreases with increasing 0, and 3) for 
low ( (< 0.003) r appears to approach a constant 
value. This latter observation is suggested in Fig- 
ure 8 by the silicone/acrylic data and by g(8) in 
equation (29) which is a constant for any given 8. 
Exploiting these observations an empirical corre- 
lation was developed using the relationship 

r=10 B ( A +g{0) = M>8) (34) 

where 


t,v , 2 \ 

’■=*?(! + '») 


u( p A 1 ' 2 


\R<ra 2 cos0/ 


1 — sin 8 


COS 0 


-1.2 a 2 + 2.2a + 0.28 

(35) 

B = 3.94 - 3.32 

(36) 

g{6) = 0.01a 2 

(37) 


A, g(0), and the function al +a 2 ) in r were 
determined empirically in terms of a. Equation 
(34) solved explicitly for the settling time t„ pro- 
duces 

t, = ^(l + a 2 )" 1 (l0V‘ + 0.01a 2 ) (38) 

The correlation of the data when represented em- 
pirically by <f> and t is illustrated in Figure 9 
with the aqueous glycerine series and data from 
other sources added [12,13]. The data of Kaukler 
agree well with the correlated data of this study 
as do those of Ross and Pline. Both researchers 
used different-sized cylinders, the latter with a 
radius three times that of the majority of this 
study. Additional data points were obtained in 
this study using distilled water and methanol in 
acrylic cylinders with and without F-723. These 
are also included on the plot. 
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The correlated settling time of equation (38) 
for surface reorientation in right circular cylin- 
ders upon release in a drop tower permits ex- 
plicit calculation to within a factor of two for 
the entire range of parameters tested; namely, 
0.65 < v < bOQ(cSt), 9.52 < R < 30(mm), 760 < 
p < 1238 {kg/m 3 ), 0.0159 < a < 0.0695(JV/m), 
and 0° < 9 < 70°. For over 85% of the data the 
accuracy of equation (38) is better than ±50%. 
For the specific case of 9 = 0° the relation sim- 
plifies to 

t, = — (50C 1J * + 0.01) (39) 

V v 

and appears accurate to ±20%. 

4.3 Sources of Scatter. 

Sources of scatter include container cleanliness, 
wall pre-wetting, photographic limitations, and 
variations in the drop release mechanism which 
can cause non-axial accelerations. To the ex- 
tent possible, these effects were minimized by the 
methods described earlier. 

The data correlating least well are the aqueous 
glycerol mixtures which have large contact an- 
gles and < 7 / p values roughly 3 times those of the 
other fluids tested. For these mixtures, numer- 
ous initial surface configurations were observed 
with only slight changes in fill technique, due to 
hysteresis effects common to large contact angle 
systems. As recognized in section 2.2, the normal 
gravity, initial shape of the liquid surface affects 
the settling time, t $ , through A* in equation (18). 
Hence the observed variability in initial shape 
likely caused the unsatisfactory correlation. 

The low viscosity fluids also produced slightly 
variable initial surface configurations depending 
on the disturbances to the cylinders as they were 
handled just prior to the drop. These fluids were 
more likely to pre-wet walls and so change the 
initial surface curvature. However, through care 
and repeated tests, this problem was virtually 
eliminated. 

4.4 Settling Time Dependence on 9. 

The general effect of varying 9 on t s may be 
seen by plotting t 9 of equation (38) against 9 . 
This is presented in Figure 10 using the proper- 
ties of silicone oil with R = 9.52mm. The curves 
of constant i/(p/R<r ) 1 ^ 2 (= Cs=o°) demonstrate 


the contact angle effects; namely, a maximum 
time required for surfaces to settle occurs at some 
intermediate value of 9 depending on the damp- 
ing of the system. The relationship is also repre- 
sented in Figure 11 where t 9 is plotted against 
v {p/ Rcr) 1 ^ 2 for cases of constant 9. It is ap- 
parent from both figures that in the vicinity of 
5 = 0° the settling times increase as 9 increases. 
This is due to the significant reduction in viscous 
damping at the wall with only marginal changes 
in interfacial potential energy. The low viscos- 
ity, large contact angle fluids ( 9 > 30°) rapidly 
established a pinned contact line and oscillations 
about this line proceeded until damped. For the 
same fluids with small contact angle, the first os- 
cillation was severe enough to cause liquid to wet 
the container wall above the final equilibrium po- 
sition. Viscous damping at the wall was accentu- 
ated in these cases due to the increased wall area 
in contact with the fluid. 

It is also apparent that t s approaches 0 as 9 
approaches 90° since the low and normal grav- 
ity surface shapes are nearly equivalent (flat). 
However, in these cases viscous damping is also 
small and wall damping is essentially nonexistent. 
Should disturbances be present (axial or non- 
axial depending on release technique, asymmetric 
pre-wetting, non-level cylinder, etc.) larger-than- 
expected settling times may result. The cases of 
9 ~ 90° are also singular points for the corre- 
lation (£ — * 00 ) and values in the near vicin- 
ity (~ 85°) may serve as a conservative estimate 
given the fact that disturbances will occur upon 
release. 

5, Concluding Remarks 

A correlation to determine the settling time 
of fluid surfaces in cylinders has been devel- 
oped for use in the design of drop tower exper- 
iments. The result also provides some guidance 
for rapidly estimating settling times given a gen- 
eral reduction in body force. This could be use- 
ful in actual space applications (e.g. when en^ _ 
gine thrust is terminated and liquid in propellant 
tanks/containers rapidly enters a low gravity en- 
vironment). Although a wide range of viscosities 
and contact angles were tested and well corre- 
lated, some caution should be taken when extrap- 
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dating the results to systems with large cylinder 
radii since this parameter was varied only over a 
factor-of- three range. 
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Table 1. Relevant fluid/interfacial properties. 


Liquid 

v± 1% 
( cSt) 

p± 5 
(kg/m 3 ) 

cr ± 5% 
(N/m) 

C 

Silicone Oil 

0.65 

760 

0.0159 

0.00146 

>? 

1.0 

816 

0.0174 

0.00222 

55 

1.5 

851 

0.0180 

0.00334 

55 

2.0 

872 

0.0187 

0.00442 

55 

5.0 

913 

0.0197 

0.0110 

55 


935 

0.0201 

0.0221 

55 


949 

0.0206 

0.0440 

55 





55 


964 



55 


967 



55 





Eth/H 2 0 20/80(% vol. ) 

1.74 

958 

0.0443 

0.00262 

” 30/70 

2.36 

937 

0.0364 

0.00388 

” 40/60 

2.64 

916 

0.0328 

0.00452 

” 50/50 

2.75 

896 

0.0308 

0.00481 

" 60/40 

2.68 

875 

0.0295 

0.00473 

” 70/30 

2.47 

854 

0.0275 

0.00446 

” 80/20 

2.36 

833 

0.0260 

0.00433 

” 100/0 

1.43 

791 

0.0227 

0.00274 

Gly/H 2 0 4.9/95.1 (% mass) 

1.04 

1010 

0.0655 

0.00132 

” 32.6/67.4 

2.25 

1072 

0.0695 

0.00286 

” 54.4/45.6 

5.36 

1127 

0.0675 

0.00710 

” 64.8/35.2 

10.6 

1155 

0.0680 

0.0142 

” 72.5/27.8 

18.5 


0.0692 

0.0247 

” 82.1/17.9 

44.4 


0.0680 

0.0606 

” 92.9/7.1 

205.0 


0.0670 

0.285 

Methanol 

0.75 

800 

0.0226 

0.00145 

Distilled Water 

1.08 

998 

0.065 

0.00137 

i-C 3 H 7 OH at 34°Ct 

1.91 

KM 

0.0200 

0.00336 

” at 42°ct 

1.15 


0.0190 

0.00207 


*, £ with 0 = 0° and R = 9.525 mm unless otherwise specified 
f, Data for Kaulkler [12] with R = 12.5mm 
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Table 2a. Results of experimentation. 


Liquid/Solid 
with Air 

Cyl. Radius ] 
(m) I 

6± 2° 
(degrees) 

Nodes 

w ± 0.2 
(eye l sec) 

t, ± 10% 
(aec) 

Si. Oil 0.65 cSt /Acrylic 

0.00952 

0° 

3 

3.1 

1.58^ 

" 1.0 

0.00952 

0° 

3 

3.0 

1.15 

” 1.5 

0.00952 

0° 

3 

3.4 

0.94 

” 2.0 

0.00952 

0° 

3 

3.7 

1.06 

” 5.0 

0.00952 

0° 

4 

4.6 

1.40 

” 10.0 

0.00952 

0° 

3 

5.3 

1.75 

” 20.0 

0.00952 

0° 

1 

5.8 

2.0 

” 50.0 

0.00952 

0° 

- 

- 

2.5 + 

” 100 

0.00952 

0° 

- 

- 

3.25+ 

” 10.0' 

0.030 

~ 0° 

3 

0.8 

5.10 

Si. Oil 0.65c5t /F-723 

0.00952 

32° 

6 

5.7 

3.38+ 

” 1.0 

0.00952 

44° 

12 

6.4 

2.75+ 

” 1.5 

0.00952 

i 46° 

8 

6.7 

1.94 

” 2.0 

0.00952 

| 51° 

6 

7.0 

1.62 

” 5.0 

0.00952 

; 59° 

9 

7.0 

1.12 

” 10.0 

0.00952 

60° 

4 

7.1 

0.66 

” 20.0 

0.00952 

66° 

4 

7.3 

0.55 

” 50.0 

0.00952 

67° 

- 

- 

0.33 

” 100 

0.00952 

68° 

- 


0.25 

” 200 

0.00952 

66° 

- 

- 

0.10 

” 500 

0.00952 

65° 

- 

- 

0.28 

” 100 

0.01904 

61° 

- 

- 

1.55 

” 200 

0.01904 

j 65° 

- 

- 

1.67 

” 500 

0.01904 

74° 

- 

- 

2.5+ 

” 10.0’ 

0.030 

I 57° 

6 

1.4 

4.72 


-f, Linearly extrapolated value 
*, Data from Ross and Pline [13] 
No measurable oscillations 


11 




Table 2b. Results of experimentation (continued). 


Liquid/Solid 

Cyl. Radius 

6± 2° 

Nodes 

u/±0.2 

t, ± 10% 

with air 

(m) 

(degrees) 


(eye/ sec) 

(sec) 

Eth-H 2 0 20/80 /Acrylic 

0.00952 

63° 

7 

10.8 

0.84 

” 30/70 

0.00952 

51° 

10 

10.0 

1.15 

” 40/60 

0.00952 

40° 

9 

7.8 

1.54 

” 50/50 

0.00952 

35° 

9 

7.8 

1.42 

” 60/40 

0.00952 

20° 

8 

6.4 

1.38 

” 70/30 

0.00952 

15° 

10 

6.1 

1.79 

” 80/20 

0.00952 

9° 

5 

5.5 

1.50 

" 100/0 

0.00952 

< 5° 

3 

4.4 

1.06 

Gly/H 2 0 4.9/95.1 /Acrylic 

0.00952 

62° 

18 

14.1 

1.30 

” 32.6/67.4 

0.00952 

66° 

22 

13.5 

1.87 

” 54.4/45.6 

0.00952 

68° 

7 

12.4 

0.75 

” 64.8/35.2 

0.00952 

68° 

8 

11.9 

0.75 

” 72.5/27.8 

0.00952 

70° 

6 

11.6 

0.58 

” 82.1/17.9 

0.00952 

57° 

- 

- 

0.23 

” 92.9/7.1 

0.00952 

57° 

- 

- 

0.12 

” 92.9/7.1 

0.01904 

57° 

- 

- 

1.50 

H 2 0/ Acrylic (Eth. rinse) 

0.00952 

67° 

11 

14.4 

0.82 

” (Gly. rinse) 

0.00952 

54° 

14 

14.0 

1.12 

1 

Methanol/ Acrylic 

0.00952 

0° 

7 

4.5 

1.53 

Methanol/F-723 

0.00952 

50° 

17 

7.9 

2.50+ 

i-C 3 H 7 OH/Glass (34°C) # 

0.0125 

0° 

2 

1.79 

14 ! 

” (42°C) # 

0.0125 

0° 

2 

1.6 

1.60 j 


+, Linearly extrapolated 
Data from Kaukler [12] 
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FIGURE 1. - NORMAL AND LOW-GRAVITY IN TERFACrAL SURFACES 
WITH AXIAL LENGTH SCALE Z IN RELATION TO R AND 0. 
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FIGURE 2. - SCHEMA! I C Or EXPERIMENTAL APPARATUS. 



Time (seconds) 

FIGURE 3. - SURFACE SETTLING HISTORIES FOR SILICONE OIL/ACRYLIC PAIR 
(0 = 0° AND R - 9.52 mm). 


13 


ORIGINAL PAGE IS 
OF POOR QUALITY 






Time (seconds) 

FIGURE *». - SELECTED SETTLING HISTORIES TO ILLUSTRATE CONTACT ANGLE EFFECTS 
FOR SILICONE OILS (R = 9.52 mm). 



Time (seconds) 

FIGURE 5. - SETTLING HISTORY FOR REFRACTIVE INDEX MATCHED FLUID 
(R - 9.52 mm). 
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FIGURE 6. - FREQUENCY OT OSCILLATION VERSUS DAMPING RATIO 
(R - 9.5? mm). 


FIGURE 7. - SETTLING TIME VERSUS DAMPING RATIO (R - 9.52 mm). 
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FIGURE 8. - N0NDIMENSI0NALI7ED SETTLING TIME VERSUS DAMPING 
RATIO (R = 9.52 mm). 
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FIGURE 9. - NONDIMENSIONAL IZED SETTLING TIME VERSUS CORRELATION FUN- 
CTION <D<r,e> FOR ENTIRE DATA SET. 



FIGURE 10. - GENERAL DEPENDENCE OF l s ON 0 AS DETERMINED BY EQUATION (38) 
FOR CASES OF CONSTANT V(p/RO) 1/? (PROPERTIES OF SILICONE OIL FOR THE 
STATED KINEMATIC VISCOSITY ARE USED WITH R = 9.52 mm). 
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FIGURE 11. - GENERAL 1)1 PI N!)[ NCI OF l $ ON f 0 = „0 AS DETER- 
MINED BY EQUATION (38) (SILICONE OIL PROPERTIES USED WITH 
R 1 O.S? mm). 
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